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ABSTRACT 


The heaving and pitching motions of ALVIN and her support catamaran are 
predicted theoretically and experimentally determined at zero speed by model 
tests. Comparison is made between theory predictions and experimental results 
for each of the vessels and then the model test results for the two vessels, 
with ALVIN in the recovery position are compared. 


The theoretical results were computed by a computer program based on the 
Korvin-Kroukovsky linear theory of ship motions in conjunction with Grim's 
added mass and damping coefficients. The experimental work with 1/20 scale 
models was performed at the M.I.T. Ship Model Towing Tank. 


It is concluded that theory predictions, correlates reasonably well with 
model test results for the catamaran. The comparisons for ALVIN were invali- 
dated by questionable experimental data and possible inapplicability of theory. 
From experimental tests, it was concluded, that the motions of the two ships 
are synchronous at wavelengths about equal to the length of the catamaran, 
that the catamaran significantly damps the heave of ALVIN, and that the ALVIN 
has little effect on.the catamaran, in the recovery position. Recommendations 
for further research in this important part of a "Deep Submergence Vehicle 
System" are presented. 
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CHAPTER I 


INTRODUCTION 


The deep submergence vehicle is but a component of a larger complex, known 
as the "Deep Submergence Vehicle System." Next to the submarine, the surface 
support unit is the most important component of this system. The support unit 
provides transportation, maintenance, launching and recovery, surface tracking, 
navigation, and communications with the submarine. 

The launch and recovery of the submarine in the open ocean presents a 
challenging problem to the designer of such a system. The essence of the prob- 
lem is how to handle the heavy and unwieldy load of the submarine safely and 
without damage in a seaway. 

Wood Hole Oceanographic Institute designed and built a catamaran support 
vessel (Figure I) during the construction of ALVIN, their deep submergence 
vehicle (Figure II). The catamaran was equipped with a cradle that could be 
lowered and raised between the hulls for launch and recovery of ALVIN. This 
was thought to be the best solution to handling of ALVIN (Dry Weight-13.6 long 
tons) in a seaway. 

During the design of the catamaran several studies were conducted on 
methods to absorb the shock expected on initial contact between the submarine 
and catamaran cradle. Because initial trials on launching and recovering a 
dummy ALVIN in a seaway, revealed that the expected shock did not exist, these 
studies were not implemented. 

Ihe catamaran concept has proven itself as an extremely successful support 


vessel in the open ocean. It was able to launch and recover ALVIN in up to 





(15) 


state 4 seas and to be a stable and seaworthy vessel in 25-ft seas 

Two main objectives comprise the scope of this investigation. The first 
is an evaluation of the accuracy of linearized theory of heaving and pitching 
motions of ALVIN and the catamaran by comparing theoretical computations to 
experimental model motion data under similar conditions. The second was to 
compare the ALVIN model motion data with the catamaran model motion data in the 
recovery position. 

The conclusions from this investigation may prove useful in the future 


design of a "Deep Submergence Vehicle System." 


FIGURE I 
CATAMARAN 
DSRVT-1 
LOA- 96 FEET 


BEAM-48 FEET 
DISPLACEMENT-350 TONS 
DRAFT, KEEL— 6 FEET 
DRAFT, PROP.—8 FT. 6 IN. 
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CATAMARAN OUTBOARD PROFILE 


(APPENDAGES WERE NOT MODELED) 


(Fig. 16 in "ALVIN, 6000 FT. Submergence Research Vehicle") 
by Mavor, etc. 
SNAME, 1966 
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ALVIN HULL LINES 
APPENDAGES WERE NOT MODELED, BUT PROPELLOR 
SHROUD AND EXTERNAL BOUYANCY BLOCKS (NOT 
by J.N. Mavor, Froehlich, Marquet, Rainnie 
SNAME, 1966 
232 


SHOWN) WERE MODELED. 
(Fig. 8 in "ALVIN, 6000 FT. Submergence Research Vehicle") 





CHAPTER II 


PROCEDURE 


General 
This study was conducted by two methods: 
1) the theoretical computations of heaving and pitching motions 
of ALVIN and of the catamaran. 
2) the experimental investigation of heaving and pitching motions 
of ALVIN and catamaran models in the M.I.T. Ship Model Towing 
Tank. 
Presented in this section are the various assumptions and procedures 
utilized in the theoretical computations and the tests conducted with the 


models. 


Theoretical 


The theoretical computations were based upon linear strip theory 


(3) 


developed by Korvin-Kroukovsky 


mass and enr UN This theory has been utilized in a digital computer 


in conjunction with Grim's results on added 


program developed at the M.I.T., Department of Naval Architecture and Marine 
Engineering by Haslum and Πτι... The analytical details of this 
linear theory as utilized in the computer program are summarized in Appendix 
B from reference (7). Use of this computer program, described in Appendix 
C, was made for the theoretical computations. 

The basic input information required for ships to the computer program 


is the; 


1) number of stations, length, and displacement 





2) 


3) 


maximum beam, sectional area coefficient, and draft for each 
station 
radius of gyration, and longitudinal location of center of 


gravity 


The assumptions made in order to submit ALVIN and catamaran input data were; 


1) 


2) 


3) 


4) 


catamaran represented as a single hull with half the displace- 
ment of the catamaran. Equivalent to assuming no hull inter- 
ference. 

the maximum beam for both ships occurred at the waterline for 
both ships 

ALVIN's shrouded propellor was represented as a solid disk 
submerged stations immediately forward of the shroud were 
entered as if a vertically solid up to the waterline from the 


maximum beam. 


The input information for waves to the computer program was: 


1) 
2) 
3) 


4) 


5) 


wave amplitude -0.125 ft., 2.5 ft. full scale 

shortest wavelength -.75 ft., 15 ft. full scale 

longest wavelength - 10 ft., 200 ft. full scale 

increment steps from shortest to longest wavelength -0.25 ft., 
5 ft. full scale 


speed of ship -0.0 


The computer output was, for each wavelength to shiplength ratio, the 


frequency of encounter, the heave amplitude (ft.), the heave amplitude non- 


dimensionalized with wave amplitude, the pitch (degrees), the pitch (radians) 


non-dimensionalized with the maximum wave slope, 27/AÀ ho. In addition, 





although not used in this investigation, the relative bow motion, the cosine 
component, the sine component, and phase angle were outputs. 
This investigation was conducted for zero forward speed of ALVIN and the 


catamaran corresponding to the launch and recovery conditions of the system. 


Experimental 


The experimental investigation was conducted in the M.I.T. Ship Model 
Towing Tank equipped with a hydraulic wave generating system. 

The tank is rectangular in cross section, 108 ft. long, 8 ft. 7 in. wide 
with a normal water depth of 4 ft. An electrically driven carriage is used 
for a point of model attachment and has an instrument patch panel for connec- 
tion to the recording equipment. A beach at the opposite end from the wave 
generating equipment essentially eliminates wave reflection. 

The waves are generated by the motion of a paddle at one end of the tank, 
which rotates about an axis at the tank bottom. The paddle fills the full 
width and depth of the tank. The paddle is powered by a high pressure 
hydraulic ram and the motion is controlled by an electronically controlled 
servo valve. For generating a train of regular waves of a given length, a 
Sinusoidal signal from a previously prepared magnetic tape is the input to 
the wavemaker control system where the desired wave amplitude is set by a 
potentiometer reducing the 5 volt peak-to-peak sinusoidal signal. A more 
complete description given by Pearlman may be found in Reference (4). 

The model dimensions were a compromise governed by the following require- 
ments: 

1) To allow the testing together of ALVIN and catamaran, the model 


scales had to be equal 





2) It was desirable to make the catamaran model less than 5 ft. 


transportation and handling 


for 


3) The ALVIN model should be as large as possible to accommodate 


the weight of instrumentation 


The scale chosen was 1/20 the full scale size. 


The full scale and model particulars are listed in Table I. 


Z > DP > Z > 8 > о 


TABLE I 


Full Scale and Model Particulars 


Length, L,ft. 
Beam(overall),ft. 
Beam(each hull), ft. 
Dratt. T,Tt. 
Displacement, tons 
Longl. C.G. aft. 

OL P-P., ft. 

Longl. radius of 


gyration, ft. 


Length, L, ft. 


Beam ,ft. 
Draft, T,ft. 


Displacement, tons 


Longl. C.G. aft. 
ο ο ο. τι. 


Longl. radius of 
gyration, ft. 


352.5 


48.0 


22.0 
8.25 
7.07 

13.6 (SW) 





4.8 
4.8 
0.7 
0.3 
0.0436 (FW) 


2.4 


1.1 
0.4125 
0.354 


0.001658 (FW) 
0.00215 (FW) 
due to approx. 
in actual model 


0.469 


0.242 





The instrumentation used to measure wave height, pitching motion ampli- 


tude, and heaving.motion amplitude is now described. 


Wave Height 


During the author's use of the tow tank, the sonic wave probe, which is 
the usual means of measuring wave height, was out of commission. The wave 
height was measured by a two wire resistance probe that was extended to half 
depth in calm water. The two wires, spaced approximately 1 in. apart with 
the water completing the circuit, is one half of a resistance bridge with the 
other half of the bridge in the carrier preamplifier. The varying water 
height, as a wave passes the probe, changes the resistance of the bridge 
providing an electric signal proportional to the change in water level. The 
electrical signal is the input to the oscillograph recorder, thus providing 
a continuous record of the wave system. The probe was linear for +1.5 in. of 
water level change and, with a scale of 10 mm/in, the oscillograph recorder 


had an accuracy of less than +0.25 mm or within 2.52. 


Pitching Motion Amplitude 


The pitch angle record was obtained using a pitch bearing mounted in 
the pitch bearing block which is rigidly attached to the model. A rotary 
variable differential transformer coil is attached to the bearing block. The 
core is attached to and rotates with the bearing. Thus the angular motion of 
the model about the point of attachement generates a signal in the coil which 
is proportional to the pitch angle and is another input to the oscillograph 
recorder. This provides a continuous record of the pitching motion as a 


function of time. 





Heaving Motion Amplitude 


The vertical oscillations of the model were measured by a linear variable 
differential transformer. The core was within an aluminum rod, which was attached 
to the model. The rod passes through a ball bearing guide, which was .attached to 
the carriage, and contains a three winding coil. 

A heave rod mounted securely in the pitch bearing and passing through 
a ball bearing guide attached to the carriage, served the purpose of allow- 
ing vertical motion only, while the pitch bearing permitted pitch angular 
rotation only. For the catamaran model; the heave rod was square in cross 
section as was its guide, to restrain the model from other motions, espe- 
cially yaw. For the ALVIN model; the tubular heave rod had, at its upper end, 
a vane that rode in a teflon slot on a teflon bearing, to restrain the model 
from other motions. The linearsyn core rod is fixed to the heave rod with a 
plexiglass bracket and thus follows the vertical motions of the model. The 
linearsyn differential transformer generates a signal proportional to the 
vertical motion of the model and, as another input to the oscillograph record- 
er, provides a continuous record of the heaving motion. 

Figure III is a picture of the catamaran model showing the instrumenta- 
tion and Figure IV is a picture of the ALVIN model showing the instrumenta- 
tion. Table II lists the instrumentation weights. FigureXXV is a typical 
oscillograph recording tape that was used for data collection. 

The models were ballasted to the correct draft with the instrument 
weights. Then the ballast was shifted to adjust for the longitudinal radius 
of gyratibn. The procedure used to determine the model radius of gyration is 


described in Appendix D. 





FIGURE III 
CATAMARAN MODEL INSTRUMENTATION 


E nd Т 





FIGURE IV 
ALVIN MODEL INSTRUMENTATION 
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TABLE II 


Instrumentation Weights 


Model Item Weight (lbs) 
Catamaran Heave Rod (Solid) 


Pitch bearing block, bearing, 
and mounting bolts 


linearsyn core rod and 
bracket 


Heave Rod (Hollow) 


Pitch heaving block, bearing, 
and mounting bolts 


linearsyn core rod and 
bracket 





Test Procedure 

To determine the heave and pitch response of ALVIN and the catamaran, 
they were first tested separately and then both models together with ALVIN 
in the recovery position between the catamaran hulls. The runs that were 
conducted are listed in Table III. 

All runs were conducted with the carriage 30 ft. away from the wavemaker 
paddle to ensure regular waves. With the calibrated wave height probe, on 
the carriage, and no models in the water, the wave amplitude settings on the 


potentiometer were determined for each wavelength. 


eque 





The catamaran model was tested for two different wave amplitudes at each 
wavelength. A single catamaran hull was then tested in the same manner. 
Because of the bow and stern symmetry of the catamaran, all runs were for 
both directly ahead and astern seas. The ALVIN was tested in directly ahead 
and directly astern seas. 

For testing together, both models were attached to the carriage with 
the heave rod and pitch bearing. This allowed both models to move freely in 
vertical motion and in pitch angular rotation. The first set of runs were 
with heave and pitch instrumentation on ALVIN. The second set of runs were 
with the heave instrumentation on the catamaran while the pitch instrumenta- 


tion on ALVIN. This method of attachment is shown in Figures V and VI. 


TABLE III 


Experimental Runs 


Model Wave Amplitude (ft.) Wavelength (ft.) 
0.125 3-0; 420,94 .5,565.05=10.0 
Catamaran 0.0625 3.0, 4.0, 4.5, 6 OF 10.0 
(Single hull) 0.125 
u. ahead 0.0625 1.0, 2.0, 3.0, 5.0, 10.0 
seas 
ay astern 0.0625 1.0, 2.0, 3.0, 5.0, 10.0 
seas 
ALVIN, Recovery 
Position, ahead 0.0625 0, 2.0, 3.0; S209) 10.0 
seas 
ALVIN, Recovery 
Position, ahead 0.0625 Sree Oe aoe 6.0, 10.0 
seas, pitching 
only 


Catamaran, Recovery 


heaving only m 3.0, 4.0, 4.5, 6.0, 10.0 


— -. 








FIGURE V 
RECOVERY POSITION MODEL INSTRUMENTATION 





FIGURE VI 
RECOVERY POSITION MODEL IN WAVES 
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CHAPTER III 


ANALYSIS 


The theoretical computed results for the catamaran were for a single hull 
only. There is no known way for the Korvin-Krovskovsky linear strip theory to 
take into account two hulls as in a catamaran configuration. Either the assump- 
tion must be made that no interference exists between the two hulls or a correc- 
tion must be made. 

The damping of heaving and pitching motion is caused by the oscillating 

body generating a traveling wave system in which energy is dissipated (2). 
With two hulls in close proximity to each other and both oscillating, there 
must be interference between the hulls. Therefore it would be erroneous to 
make the assumption of no interference and accordingly a way to correct the 
theoretical results was sought. 

Because there was no theory available, empirical methods were the only 
means available and the procedure used to correct the theoretical results is 
described herein. 

The non-dimensional experimental heave and pitch motion amplitudes for the 
catamaran and a single hull were plotted for two different wave heights, 

(Figures VII, VIII, IX, and X.) From the data obtained for the two different wave 
heights, at each wavelength, a mean value was obtained by interpolation for 
both the case of the single hull and the catamaran. These points are plotted 
in Figures XI and XII. The ratio of the catamaran to the single hull was 
computed for the different wavelengths and applied as a correction to the 
theoretical non-dimensional motion amplitudes for the single hull. (See 


Table IV). 
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TABLE IV 


Correction for Single Hull to Catamaran 


Mean | Mean Cat. 
A(f£t)} Single] Cat. Single | Theory 


(heave amplitude) 













Corrected 
Amplitudes 






©1200} +.1200 





(pitch amplitude) 


.0501| 1.67 „0749 

.1600| 2.42 .1300 
1.87 . 3080 
1.31 „3525 
1.32 „4760 
1.25 «6325 
1.228 . 7480 
1. 232 8430 
1.242 «9200 
1.262 





. 9850 
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CHAPTER IV 


RESULTS 


The results are presented in the form of non-dimensional motion amplitudes 
versus wavelength and wavelength to shiplength ratio, Figures XIII-XXIV. Heave, 
Zo? is non-dimensionalized with the wave height, ho» and pitch, са in radians, 
by the maximum waveslope, 2T/À ho: Full lines represent theoretically computed 
motion amplitudes and experimentally determined motion amplitudes are represented 
by circled points. 

Figures XIII and XIV compare catamaran theory predictions with experimental 
results. Figures XV and XVI show ALVIN comparisons in Heave and Figures XVII and 
XVIII show ALVIN comparisons in Pitch. Figure XIX is a comparison of ALVIN with 
and without heave instrumentation. Figures XX and XXI are comparisons of the 
two vessels in the recovery position. Figures XXII and XXIII are comparisons of 
ALVIN alone and in recovery position. Figure XXIV compares the catamaran alone 


and with ALVIN in recovery position for heave motion. 
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CHAPTER V 


DISCUSSION OF RESULTS 


For the catamaran, after the theoretically computed heaving motion 
amplitudes had been corrected, the plotted curve correctly predicted the trend 
of experimental results. While theory ТЕН somewhat higher at shorter 
and longer wavelengths, in the range of 0.9 to 1.5 wavelength to shiplength 
ratio there was reasonably close agreement with experimental results. (Figure 
XIII). 

The curve of corrected theoretical pitching amplitudes correctly predicted 
the trend but was about 20% higher than experimentally determined results. 
(Figure XIV). 

Instrumenting the small ALVIN model was extremely difficult because the 
weight of instrumentation was of the order of magnitude of weight of the model. 
A hollow heave rod of tubular aluminum was substituted for the solid square 
steel heave rod. However, during the model testing for motion Er. the 
light weight of the instrumented ALVIN was not sufficient to overcome the small 
friction in the guides of the heave rod and linearsyn core rod. For many runs 
no vertical oscillation was recorded at all. The runs where heave motion was 
recorded showed intermediate binding so that all experimental heave data of 
ALVIN was considered questionable. 

Because of the undesired restraint on vertical motion, the model responded 
to the wave action with amplified pitch motion. An example of this is shown 


in Figure XIX, which compares the pitching motion with and without the linearsyn 


зара 





rod. This proves that pitching motion was artificially increased due to this 
unwanted restraint. Another piece of evidence is the fact that at 5.0 ft. wave- 
lengths and 10.0 ft. wavelengths pitch exceeded by a factor of 2-5 the maximum 
waveslope. At pitch angles about equal to 5 degrees bow up the center of 
buoyancy shift was great enough to cause an instability in the model and it 
would flip stern down. For this reason, the maximum wave amplitude that could 
be used with ALVIN was 0.75 in. crest to trough. 

Comparisons with theoretically computed heave and pitch motions amplitudes 
were nevertheless made, albeit somewhat questionable. Theoretical results were 
higher for both directly ahead and astern seas. Also the curve of theoretical 
results of heaving amplitudes predicted a resonance at the 2.25 ft. wavelength 
that was not verified by the experimental data. (Figures XV and XVI). For 
ALVIN pitching motions the experimentally determined results crossed the 
theoretical curve at about 4.0 ft. wavelength in both directly ahead and astern 
seas where the restraint in vertical motion caused greater pitching motion of 
the model. (Figures XVII and XVIII). 

The recovery position motion data shows that for heaving amplitudes the 
two vessels first follow the same trend with ALVIN heaving slightly more than 
the catamaran at wavelengths shorter than the catamaran length then about the 
Same at a wavelength 90% of the catamaran length. Then the curves of heaving 
amplitude diverge with the catamaran heaving amplitude greater. (Figure XX). 
This is contrary to the author's observations during the runs, when it was 
directly observed at the 4.5 ft. wavelength the ships seemed to move in perfect 
Synchronism and exhibited similar, though less perfect, synchronism as the wave- 


length increased. The discrepancy between measured results and observation may 
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be caused by binding of the heave instrumentation on the ALVIN model. 

For the recovery position pitching motion, the rapid divergence of ALVIN 
at longer wavelengths is felt to be, again, due to the restraint in heave. 
(Figure XXI). 

ALVIN's effect on catamaran heaving motion appears to be slight and 
highly frequency (wavelength) dependent. (Figure XXII). 

ALVIN's heaving motion amplitude appears to be significantly damped 
(Figure XXIII). 

ALVIN's pitching motion, however, is amplified at the shorter wavelengths 
with a resonance at a wavelength 62% of the catamaran, and damped slightly 
at a wavelength 90% of the catamaran length. (Figure XIV). This is, again, 
inconclusive because of the ALVIN model instrumentation. 

Reasons for Discrepancies 

The experimental methods and the current status of theoretical knowledge 
possess many sources of errors. To provide a partial explanation of the 
differences in this investigation between theory and experiment and the 
differences between the models, these discrepancies are discussed. 

For the experimental work, the following sources of discrepancies exist: 

a) Wave height measurement is one of the greatest sources of errors 


(7). 


in tow tank work The potentiometer settings for wave- 
heights were determined with no models in the tank to 
eliminate the model effect on the wave height measured. How- 


ever, as the potentiometer settings were quite delicate, this 


may have led to errors. 


Loc 





b) 


с) 


d) 


e) 


f) 


Wall effects, or reflection from the walls back onto the model 
of waves generated by the model motion at zero forward speed 
and at small wavelengths significantly affected the ee 
Undesired instrument effects on the model may have affected the 
results, as was previously discussed for the ALVIN model. 
Instrument accuracy was about +2.54. 

Extraction of recorded data on oscillograph recorder tapes may 
have introduced errors. 

During the 1.0 ft. wavelength generation, the hydraulic ram 
pushing the paddle was observed to loosen on its foundation 
and a modulated wave record resulted. This may have affected 


other wavelengths, but oscillograph records appeared to be 


regular. 


For the theoretical computations, the following assumptions limited the 


applicability of the theory to ALVIN and the catamaran: 


a) 


b) 


linearized theory assumes ship section lines vertical at the 

waterline. Neither of the models have vertical sides, espe- 

cially in the case of ALVIN. 

the computation of damping in linear theory is by accounting 

for the energy dissipated in gravity waves generated when the 

ship oscillates. Frictional, separation, and circulation 

phenomena are assumed to be negligible, which is generally valid 
(5) 


for surface ships o This assumption is erroneous for ALVIN 


where the above phenomena are a considerable influence. 


ag 





c) 


d) 


e) 


linear theory assumes small motions and that coefficients of 
motion are constant throughout the whole cycle of motion. 

This is not true for larger motions. The restoring force 
coefficients are correct only for vertical sides at the water- 
line. For the ALVIN model with a decreasing and increasing 
waterplane area as the model oscillates up and down, this is 
certainly not true. 

end effects or three-dimensional effects are responsible for 
large discrepancies when large section areas are present at the 
ends of the model. The sectional area curve must go "grace- 
fully" to zero at the bow and stern in order for the force 
contributions of these sections to zero. This does not affect 
the heave coefficients seriously but the ends of the ship have 
a double effect on pitching and therefore significant errors 
may result in pitching Da) This is serious inadequacy 
in applicability of the theory to ALVIN. 

linear strip theory has a basic assumption that beam and draft 
are small compared to length and uses a basically two-dimen- 
sional T This is another serious inadequacy in the 


applicability of the theory to ALVIN, which has a length to 


beam ratio of 2.6 and a length to draft ratio of 3.1. 
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CHAPTER VI 


CONCLUSIONS 


The conclusions from the preceding chapters are summarized. 


1, 


The catamaran corrected theoretical results, in spite of discrep- 
ancies noted previously, predicts and agrees reasonably well with 
the experimental results. 

The Korvin-Krovkorsky linear theory, with Grim's added mass and 
damping, applicability to ALVIN is questionable but this conclusion 
cannot be substantiated because of inaccurate experimental data. 
Because of the inapplicability of theory and the questionable 
experimental data, a valid comparison could not be made between 
ALVIN theoretical predictions and experimental results. 

ALVIN and the catamaran respond in synchronous motion at wavelengths 
about equal to the length of the catamaran, but this conclusion is 
questionable because of ALVIN experimental data. 

ALVIN heaving motion is appreciably damped in the recovery position 
between the catamaran hulls, butthis conclusion is questionable 
because of ALVIN experimental data. 

The catamaran heaving motion is not significantly affected by ALVIN 


in the recovery position. 
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CHAPTER VII 


RECOMMENDATIONS 


Much research remains to be done in the design of surface support vessels 


for lauching and recovery of deep submersibles in the open ocean. Therefore it 


seems important that further theoretical and experimental work be undertaken in 


this important part of the "Deep Submergence Vehicle System." 


Recommendations arising from this investigation are that; 


1) 


2) 


3) 


4) 


the theory be modified to extend its applicability to deep sub- 
mersibles in surfaced condition similar to the hull form of ALVIN. 


There exists in the literature the more general analysis of πο ο 


(18) 


and Bermejo that accounts for section shape in more detail. 

the comparisons attempted in this investigation be repeated with 
special attention given to devising better heave instrumentation 
for ALVIN or larger scale models be used. 

this investigation should be extended by experimentally determining 
phase angles of motion and comparing them with the theoretical 
computed phase angles (that were not used in this investigation). 


an experimental investigation be made using different support vessel 


hull forms to determine the optimum hull form for open ocean. 


ο 
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APPENDIX A 


NOMENCLATURE 


coefficient of equation of motion 
coefficient of equation of motion 
coefficient of equation of motion 
coefficient of equation of motion 
coefficient of equation of motion 
coefficient of equation of motion 
coefficient of equation of motion 
coefficient of equation of motion 
coefficient of equation of motion 
coefficient of equation of motion 
total vertical (heaving) force 


amplitude of time-varying heaving force 


complex vertical (heaving) force 

gravitational acceleration 

coefficient of equation of motion 

coefficient of equation of motion 

amplitude of sinusoidal wave (half-wave height) 
maximum wave slope, (2π/λ ho) 


longitudinal radius of gyration 


length of ship or model 
mass of ship or model 


total pitching moment 


SS 





= | 


N(ë) 


= amplitude of time-varying pitching moment 


= complex pitching moment 


il 


sectional damping coefficient 
= time 
= speed of ship or model 


= heaving motion of center of gravity of ship or model 
= heaving velocity of center of gravity of ship or model 


= heaving acceleration of center of gravity of ship or model 


= amplitude of heaving motion 


= complex heaving motion 


Greek Letters 


theoretically computed heaving phase angle 
displacement of ship or model 

theoretically computed pitching phase angle 
pitching displacement 


pitching velocity 


pitching accelerations 


amplitude of pitching motion 


complex pitching motion 
wavelength 

water density 

phase angle of heaving force 
phase angle of pitching moment 
frequency of encounter 
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APPENDIX B 


ANALYTICAL DETAILS OF THE LINEAR THEORY OF SHIP MOTIONS 


In order to achieve completeness and simultaneously provide easy reference, 
the fundamental analytical expressions of the Korvin-Kroukovsky method, as 
utilized in the computer program, are summarized in this Appendix. Following 
the nomenclature and definitions adopted in (8), the coupled set of linear 


differential equations describing the two-degree-of-freedom ship system, takes 


the form, 
aZ + bz + cz + dO + e0 + gÔ = F exp (iu, t) 
Т | . = (1) 
АО + BO + CO + DZ + Ez + Gz = M exp (iu t) 


5 


The above equations result from equilibrium considerations of the hydro- 
dynamic forces and moments called into play by the ship's oscillations in the 
plane of symmetry, when meeting head or astern regular waves. For this reason 
the analysis ignores steady, continuously acting forces due to buoyancy, gravity 
and suction pressures. Following the principles of classical dynamics, these 
forces and moments are obtained by applying Newton's Second Law of Motion to 
both translatory and rotational displacements of the body's center of gravity. 

The wave induced excitation force and moment may be thought of as being 
imposed on a fully restrained ship and appear on the RHS of the set (1). 

They have the useful property that they are functions of the wave elevation and 
its two first time derivatives, thereby allowing ease in algebraic πος. 


They are defined in complex notation as 


| 
I 
I 


exp (іш t) F. exp (-i0) exp (iu t) F. exp [1 (w_t-0)] 


(2) 


c4 
и 


ехр (iu t) M exp (-iT) exp (io t) M. exp [i(o t-1)] 
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The differential exciting force acting on a control section distant Ë from 


the origin of the moving coordinate system (ship's C.G.), has been given in 


(8) in the simplified form, 


dF 

dx 
where, 

9] 
and, 

Ф 








аЕ dF 
= cos wt + mm sin we 
2Thc 
2 2T 
(1%, біп az + $, 5 x cos 215) exp (- = боз wit 
(3) 
2Thc 
. τπτ 21 | 
11%, сов 215 - Ó$, 5 - зіп 215) ехр (- +*)] sin ш © 
ат вс“ 
hog ВЖ - nu (ρ Sk,k,) 
A 
d(p Sk,k,) 
2 4 
NG) - V — 


while the differential exciting moment of this force about the C.G. is given by 


dF 


dx ἀξ. 


Integration of the above two quantities over the ship length results in 


the values of the total time-varying exciting force and moment, which are con- 


sidered as the real parts of the expressions (2). 


F 


1 


ШІ? 
2 Ῥ 


= F cos(w_t-0) 


к. cos t + F 
e 


2 


2 cos[w t-arctan 


Thus, 


sin wt M M. cos ü t + M 
e 1 e 


= 


Е М 
2 _ [2 2 2 
г! апд ο. cos[w t arctan m 


M соз(ш Е-т) 
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2 2 n we 


ll 
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The analysis of the forces and moments which correspond to the ship's free 
oscillations in calm water yields terms which appear on the LHS of the set (]) 
and are proportional to the instantaneous heaving and pitching displacement, 
velocities and accelerations. All twelve coefficients of the above terms are 
independent of the speed per se and the body's space orientation, but depend 
on the frequency of encounter, with the exception of c and G. The final 
expressions for the coefficients of the equation of motion used in the machine 


computations are listed below: 


a = m+tp | (Sk,k,) dE 
b = f N(E) de 

c = pg f Be dé 

d = pf (Skk,) Е d£ 


е = J N(E) € dé - 2V p f (Sk.k,) d£ - V p f à (Sk k,)/d£ E d£ 


g = pg/[B*Edt - V f N(E) ac (5) 
A= J+pf (Sk,k,) E^ de 

в = |6) E^ dE - 2ν p f (Sk,k,) E dE -VpJfd (Sk,k,)/d€ E* ἀξ 
C = pg f B*E ag - v f N(£) E at * v? p f à (Skk,)/at E de 

р = pf (sk,k,) E de 


Е = [ NE) EdE-Vp/d (Sk,k,)/dE E dE 


G = pg f B* E dE 


Assuming now that sufficient time has elapsed for any transient distur- 
bances to be damped out, we seek particular solutions of the nonhomogeneous set 


(1), which correspond to the steady-state responses of the system. Since the 
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forcing functions are sinusoidal and the system is linear and time-invariant, 
we expect that any response will also be sinusoidal of the same frequency as 
the excitation and with generally different amplitude and phase. We therefore 


assume solutions of the forn, 


z(t) = z = z exp (io t 
Ы. (6) 
апа O(t) = Ө = О exp (iw_t) 
with the arbitrary definitions, 
2 = 2 ехр(-18) 
5. (7) 
and О = 9 exp (-ie) 


Upon substitution in the original equations (1), a conscientious algebraic 
manipulation yields the following expressions for the "complex" heaving and 
pitching amplitudes, 


7 = EM апа 


PS-QR O = PS-0R (8) 
where, 
F = F (сов O -i sin 0) 
M = M (cos 1 -i sin T) 
2 А 
Р = (c-au ) tib o 
e e 
(9) 
S = ЕЕ )+iBw 
e e 
Wc (RE EN T 
е е 
2 
R = (G-Du ) + ií E ω 
e e 


It follows from (8) and (9) that the "complex'' amplitudes may be expressed 


as, 


5] 





м | 


2. 
J 25242 | E 
2 21 +z, exp[-i arctan =! 


ZJ - 12 + 
(10) 
= 2 2 9, 
and © 9, - 1 O, = 9, +0, exp[-i eg 
Considering the real parts of (6), we finally obtain, 
z = R, Ζ exp (iw, t) 
23 Ko 
= Re |z. +z expli(w t-arctan —-] (11) 
1 2 € zi 
= 20 соз(ш Е-9) 
and Ө = R, Ө exp (iu, t) 
Θ 
de 202 _2 
= Re 9, +0, ехр[і(ш е arctan Ө, 
(12) 


= 9, cos(w_t-E) 


2698 





APPENDIX C 


DESCRIPTION OF COMPUTER PROGRAM 


Table of Contents 


Nomenclature of Computer Program 
Description 

Main Program 

Subroutine ADMAB 
Subroutine COEFF 
Subroutine EXCITE 
Subroutine MOTION 
Subroutine BENDISH 
subroutine SIMPS 
Description of Input Data 
Description of Output Data 
Input and Output Data 


Listing of Program and Subroutines 


ууш 





ABAR (I) 
ADDA (1) 
ALPHA 
AVETHE 
BEEB (1) 
BETA 
BMIMAG 
BMNULL 
BMREAL 
BPL 
BSTAR (I) 
CGGC(T) 
CTFST 
CW 
CXFST 
DELTA 
DELV 
DELWL 
DISPL 
DIX(I) 
DMASS (1) 
DRAFT (1) 
DWEIGH (1) 


EMNULL 


NOMENCLATURE OF COMPUTER PROGRAMS 


see subroutine ADMAB 

see subroutine COEFF 

phase lag of shearing force 
average pitch amplitude 

see subroutine COEFF 

phase lag of bending moment 

sine term of bending moment 
bending moment amplitude 

cosine term of bending moment 
ship (model) length between perpendiculars 
full beam (breadth) of ship (model) station at DWL 
see subroutine COEFF 

see subroutine BENDSH 

wave celerity 

cosine term of exciting force 
phase lag of heaving motion 
increment of ship or model speed 
increment of wavelength 

ship (model) displacement 
d(Sk,k,)/dE 

station mass 

station depth (draft) 

station weight 


pitching moment amplitude 


ue 


ENOX(I) 
EPSIL 
FNULL 
FROUDE 
GAMMA 
GRAV 

H 

I 
INCRES 


KRIT 


MAXKRI 


MINKRI 


OMEGA 
OMEGAE 
ЕТ 


QUANT (Т) 


RO 

SECOE (1) 
SHIMAG 
SHNULL 
SHREAL 
SIGMA 


SKLAM 


station damping coefficient 
phase lag of pitching motion 
heaving force amplitude 
Froude number 

specific weight of water 
gravitational acceleration 
interval of Simpson's integration 
program variable 

see subroutine BENDSH 

see subroutine BENDSH 

no. of ship or model stations 
see subroutine BENDSH 

see subroutine BENDSH 

highest station number 
absolute wave frequency 
frequency of encounter 

т = 3.1415926 


Sk,k, 


density of water 

section area coefficient 

sine term of shearing force 
shearing force amplitude 
cosine term of shearing force 
phase lag of heaving force 


see subroutine EXCITE and BENDSH 
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STFST 
SXFST 
SYMPS 
TAU 

TI(I) 
TIMAG 
TMASS 
TNULL 
TR(I) 
TREAL 
UR(I) 


UI(I) 


VMIN 
WÀ 
WAVEN 
WL 


XI (1) 


YNERT 
ZIMAG 
ZNULL 


ZREAL 


see subroutine BENDSH 

sine term of heaving force 

area under any curve (see subroutine SIMPS) 
phase lag of pitching moment 

see subroutine MOTION 

sine term of pitching motion 

total mass of ship or model 

amplitude of pitching motion 

see subroutine MOTION 

cosine term of pitching motion 

see subroutine MOTION 

see subroutine MOTION 

ship or model speed 

maximum speed of ship or model 

minimum speed of ship or model 

wave amplitude 

wave number 

wavelength 

longitudinal coordinate of moving (ship) system 
see subroutine SIMPS 

weight moment of inertia of ship or model 
sine term of heaving motion 

amplitude of heaving motion 


cosine term of heaving motion 


Шер 
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Description 


This computer program was originally written by K. Haslum l 2598 subse- 


quently debugged and completed with the aid of L. προ, 7. The program 
has the ability to compute shearing forces and bending moments at any station 
of a ship or model encountering regular waves. Subroutine BENDSH which 
essentially performs the above task, was always bypassed in the main computa- 
tion cycle through the inclusion of a proper control card. 

This program is a computer alogorithm of the step-by-step procedure 
outlined in (8) for the calculation of ship (model) motions and bending moments 
on the basis of the Korvin-Krovkovsky linear theory of pitching and heaving. 


Since the program is extensively described in (19), only a brief description 


of its MAIN program and Subroutines is given. 
MAIN Program 


The purpose of the MAIN program is to perform certain minor calculations, 
control the computation cycle from subroutine to subroutine and process all 
input and output data. 

The program output capabilities come under three options, of which the 
present investigation has made sole use of option No. 1. The program computes 
motions only and assumes knowledge of radius of gyration. 

The main variable for the program is the wavelength and the basic computa- 
tion cycle repeats the fundamental loop, first for different wavelengths and 
then returns if the ship speed is changed. For this investigation the ship 
Speed was zero. For the complete calculation the program prints out the wave- 


length (WL), ship speed (V), and frequency of encounter (OMAGAE). The units 


t m 





are consistent and correspond to those chosen in the first READ statement. 


(See Description of Input Data). 


Subroutine ADMAB 

This subroutine, obtained from the staff of the Davidson Laboratory of 
the Stevens Institute of Technology, calculates according to Grim's ον. 

a) the added mass per unit length at each station, which is designated 

as QUANT(I), and, 

b) the ratio of the amplitude of the emitted wave to the amplitude of 

heave at each station, which is designated as ABAR(I). 

Although individually checked for normal Lewis ship sections, the range 
of the numerical applicability of the subroutine has not been thoroughly 
examined. 

Subroutine COEFF 
This subroutine evaluates numerically the expressions (5) for the coeffi- 


cients of the equations of motion. The correspondence between coefficients 


and program variables is shown below: 


анъ» с ADDA(1) ВЕЕВ(1) CGGC(1) 
deg В ADDA(2)  BEEB(2)  CGGC(2) 
DEG ADDA(3)  BEEB(3)  CGGC(3) 
A B C АррА (4) ВЕЕВ(4) CGGC(4) 


The damping coefficient at each station, denoted by ENOX(I) is obtained 


from: 


de^ e 


ω 
ee 
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where A is the quantity designated in subroutine ADMAG as ABAR(I) and the 
remaining symbols have their usual meaning. The numerical integrations are 
performed by subroutine SIMPS with an integration interval H = DX1. 
Subroutine EXCITE 

Subroutine EXCITE computes the excitation force and moment as discussed 
in Appendix B. The elemental excitation force at a given ship station is 
computed from equation (3) by splitting the above expression into two parts. 
The first designated CXFST(I) is in phase with cos wt and is given by the 
first bracketed term of equation (3). The second designates SXFST(I) is in 
phase with sin wt and is given by the second bracketed term of equation (3). 
The quantities Фі апа $^ of equation (3) are correspondingly referred to as 
FKLAM and SKLAM(I). SKLAM(I) is conveniently subscripted for subsequent 
utilization in subroutine BENDSH if the computation of bending moments is also 
required. The total excitation force and moment are finally obtained from 
equation (4) by integrating over the ship length the appropriate elemental 
values using subroutine SIMPS with H = DXI. 
Subroutine MOTION 

Whenever the main program calls this subroutine, the excitation terms and 
the coefficients of the equations of motion (1) are already available and hence 
the heaving and pitching motion amplitudes and phase angles can readily be 
determined. Subroutine MOTION performs this task by first expressing the set 


(9) in the following form: 


F = F. + iF, 
M = М. + iM. 
P = P. + iP. 
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ο. ΠΙΟ. 
R = R. + iR, 
S = S + 15. 

r 1 


The various products in the set (8) are then evaluated as shown below: 


MQ = (MQ MQ) + iM Q+ Q) 
FS = ES) + 1(F_S,+F,S_) 
QR = (QR QR) + i(Q R +Q R ) 
PS = CAE E) E iP S +P S) 
FR = (F R -F R) + i(F R +F R) 
MP = (MP _-M,P,) + і(МР,+МР) 


so that (8) can finally be computed from 


ZREAL - iZIMAG 


N 
1 


© = TREAL - iTIMAG 
where ZREAL та 
ο 
ZIMG = = E E 
Ge 
inpar > 
ci +p 
TIMAG = = Ξ = 
Стр 
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and A = Μα. - MQ, - F S. 
= dr 1 mr » t 
ы o Е ο. Е ο 
D = Q R, + Q R. - Р 5, 


Е = FR = F IR, = M P. 


Ε "ΕΕ τη». 
ICE ri 


, 
E 
rg 


The amplitudes and phase angles of heaving 


computed with equations (11) and (12), i.e., 


z 


and pitching motions are finally 


' | (ZREAL)? + (ZIMAC)? 





8 = 
ο 
ZIMAG 
and 6 arc tan {- zur) 
Ν TIMAG 
€ = arc tan Í TREAL? 


| (TREAL)2 + (TIMAG)2 


The correspondence between program variables and the quantities referred 


to above is as follows: 


D NU TR(1) 

E ΝΕ; TR(2) 

Q Q, 4 TR(3) 

κ. R, TR(4) 

M. Μι TR(5) 

5. 8 TR(6) 
-{ Ес 


=o 


ΤΙ (1). 
тт (2) 
тт (3) 
TL(4) 
TI(5) 


TI(6) 





and 


E F UR(3) UI (3) 


С р - UR(4) UI(4) 
A B |UR(5) UI(5) 


It should be noted that TR(2), TI(2), TR(5) and TI(5) are furnished by sub- 
routine EXCITE. 
Subroutine BENDSH 

For a complete description of this subroutine the interested reader is 
referred to (19), which presents the basic analytical treatment involved in 
the computation of ship bending moments in regular waves. Since this subroutine 
is at present not used by the computer program under description, we shall 
refrain from further details. 
Subroutine SIMPS 

This subroutine is an integrator and calculates numerically various 
integrals and areas under curves by using Simpson's rule. At any stage of the 
computation cycle, the number of ordinates furnished to this subroutine must be 


odd. Thus any integral is evaluated in accord with: 


b 


f f(x)dx = Е (Y, + 4Y, + 2Y 


+ о о о + 
К 3 AYN T^Y 


] 


3 N+1 


where Y(I) are the integrand ordinates and I rises from 1 to J. Y(I), H and 


J must be preset appropriately to ensure that J is an odd number. 


Description of Input Data 
The MAIN program includes six READ statements, numbered 1001 - 1006, 


whose purpose is to (a) transmit the basic input information required for the 
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computations and (b) incorporate the control cards which decide on the path of 
the computation and subsequently the form in which the output information is 
to be given. 
Unless specifically stated in the following, numerical data must be of 
the floating point type. Any consistent system of units may be used, the 
decision being made by the user for the input variables of the first READ state- 
ment (No. 1001). At the end of the FORTRAN or binary deck and after the *DATA 
card, the basic information for the computations should be given in the follow- 
ing order, one card per READ statement, unless otherwise stated: 
a) READ statement 1001, (I10,4F10.4) 
1) N - number of stations the ship is divided into. 
This number is also equal to the highest 
station number when the latter runs from 0 
at the F.P. through N at the A.P. N must be 
an integer and also, due to the structure of 
subroutine SIMPS, it must be even. Unless the 
DIMENSION statements are increased, N should 


always be less than or equal to 20. 


2) СЕРІ, = length between perpendiculars of ship or model. 
3) GAMMA = specific weight of water 

4) GRAV = gravitational acceleration 

5) DISPL = displacement of ship or model 


b) READ statement 1002, (3F10.4) 


This statement requires a number of cards equal to N + 1 giving, for 


each station starting from the F.P., the following information per card: 
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full beam (breadth) at that station 


1) BSTAR(I) 


2) SECOE(I) sectional area at that station divided by the 
area of the circumscribing rectangle at that 
station. 

3) DRAFT(I) = actual depth (draft) of section at that station. 


c) READ statement 1003, (2F20.4) 


This card must always be included since it is a control card which dictates 


the computer to bypass subroutine BENDSH. The following information is assumed 


known and must, therefore, be given: 


1) RADGYR = radius of gyration of ship or model 
2) XI(I) = distance of center of gravity of ship or model from 
the F.P. 


To compute bending moment, use a blank card for this READ statement. 
d) READ statement 1004 

This statement is never encountered in the computations described 

herein; therefore no input card should be included for this statement. 
e) READ statement 1005, (3110) 

The input card to this statement must always be included since the 
printed information serves as a controller to the program. Ordinarily, MINKRI, 
MAXCRI and INCRES are input variables for subroutine BENDSH which, however, is 
never used in the present computations. For this reason, the following 
numerical values must be given so as to instruct the computer to bypass BENDSH 


and avoid related calculations: 


MINKRI = zu 
MAXKRI = 1 
INCRES = 1 


Tor 





Е) READ statement 1006, (7F10.4) 


1) 
2) 
3) 
4) 


5) 


6) 


7) 


WA 


SWL 


BWL 


DELWL 


VMIN 


VMAx 


DELV 


il 


wave amplitude, 

shortest wave length, 

longest wave length, 

increment in wave length from SWL to BWL, 


slowest speed, 


(units consistent with READ statement 
1001). 


highest speed, 


increment in speed from VMIN to VMAX. 


Description of Output Data 


The program prints out the absolute heave and pitch amplitudes, the non- 


dimensional amplitudes, the cosine term, the sine term, and the phase angles. 


The particular modification of this program used, computes relative bow motion 


(REBM). Only the heaving and pitching non-dimensionalized amplitudes were 


utilized in this investigation. 
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Input and Output Data 


Catamaran Model - Directly Ahead Seas 


Input 
EXECUTION 
20 4.8000 62.4000 32.2000 48.8050 
0. 0. 0. 
„5200 „7625 „1875 
„6156 „7940 „2625 
„6750 «8290 «2781 
«7000 „8350 „3000 
„7000 „8350 „3000 
„7000 „8350 „3000 
„7000 „8350 „3000 
„1000 „8350 „3000 
„7000 „8350 „3000 
„7000 „8350 „3000 
„7000 „8350 „3000 
„7000 „8350 „3000 
„7000 „8350 „3000 
. T000 .8350 .3000 
.7000 .8350 . 3000 
. 7000 .8350 „3000 
„6750 „8290 „2781 
„6156 „7940 „2625 
«5200 „1625 „1875 
0. 0. 0. 
ENTER RADGYR AND CTHER DATA NOW 
.9065 2.4000 
0 0 0 
.1250 .7500 10.0000 „2500 0. 0. 
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ALVIN Model - Directly Ahead Seas 
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EXECLTICN 
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0. 0. 0. 
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.4020 .8440 .34€5 
.4015 ο8350 «3410 
„3957 «8250 „3360 
„3863 „8264 „3300 
«3765 .8228 „3225 
3632 „8210 „3165 
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Listing of Program and Subroutines 


THE DYNAMIC BENDING MOMENT IN REGULAR WAVES 


DIMENSION Y(21),TR(6),TI(6) ,ADDA(4) , BEEB(4) , CGGC(4) ,UR(6) ,UI(6), 
1DMASS(21),QUANT(21),SKLAM(21) , BSTAR(21) , CXFST(21) , SXFST(21), 
2CTFST(21),STFST(21),XI (21) , DIX(21), ENOXI (21) , DRAFT (21) ,DWEIGH(21) , 
3SECOE (21) , ABAR(21) 

COMMON Y,J,SYMPS,DXI,ADDA, BEEB , CGGC, ZREAL,ZIMAG, TREAL, TIMAG, ZNULL, 
1TNULL,DELTA,EPSIL, TR,TI,V,DMASS, QUANT, OMEGAE , SKLAM, KRIT,RO,GRAV, 
2BSTAR, CXFST, SXFST, ALPHA, SHNULL, XI , BETA, BMNULL , GAMMA ,DIX,M,WA ,WAVEN 
3,CW, ENOXI, SIGMA, TAU, FNULL, EMNULL, DRAFT , DWEIGH, SECOE, TMASS ,N,UR,UI, 
4ABAR,PI,SHREAL, SHIMAG, BMREAL, BMIMAG, YNERT , BPL 
FORMAT (110,4F10.4) 

FORMAT (3F10.4) 

FORMAT (F10.4) 

CONTINUE 

FORMAT (4F20. 3) 

FORMAT (2F20.4) 

FORMAT(7F10.4) 

FORMAT(3110) 

FORMAT (120) 

FORMAT (F9.4,2F10.4, 3F9.4,F8.4) 

FORMAT(55H FROUDE WL/BPL OMEGAE HEAVE PITCH BEND) 
FORMAT (64x18) 

SHIP-AND WATER-CHARACTERISTICS 
PI=3.1415926 

READ 900,N,BPL,GAMMA,GRAV,DISPL 

PRINT 900,N,BPL,GAMMA,GRAV ,DISPL 
RO=GAMMA/GRAV 

FN=N 

DXI=BPL/FN 

M=N+1 

READ 901, (BSTAR(I) ,SECOE(I) ,DRAFT(I) , I=I,M) 
PRINT 901, (BSTAR(I),SECOE(I),DRAFT(I),I=1,M) 
J=M 

TMASS=DISPL/GRAV 

FORMAT(32H ENTER RADGYR AND OTHER DATA NOW) 
PRINT 916 = 

READ 905, RADGYR,XI(1) 

PRINT 905, .RADGYR, XI(1) 

IF (RADGYR) .2,2,5 

DO 3 I=1,M 

READ 902, DWEIGH(I) 

PRINT 902, DWEIGH(I) 
DMASS(I)=DWEIGH (I) /GRAV 

HOMENT=0.0 

DO 4 I=1,M 

L=I-1 

FL=L 

HOMENT=HOMENT+DWEIGH (I) *FL 

XI (1) =DXI*HOMENT/DISPL 
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DO 6 I=2,M 

τι 

FL=L 

XI(I)=XI(1)-DXI#FL 

IF (RADGYR)7,7,9 

YNERT=0.0 

DO 8 I=1,M 
YNERT=YNERT+DWEIGH (1) * (XI (1) *XI(1)+(DXI*DXI(12.)) 
YNERT=DISPL*RADGYR*RADGYR 

READ 907, MINKRI,MAXKRI, INCRES 

PRINT 907, MINKRI,MAXKRI, INCRES 

READ 906, WA,SWL, BWL,DELWL, VMIN, VMAX, DELV 
PRINT 906, WA,SWL,BWL,DELWL, VMIN, VMAX, DELV 
PRINT. 909 

V=VMIN 

WL=SWL 

IF(MINKRI) 14,14,13 

KRIT=MINKRI 

WAVEN=2.*PI/WL 

CW=SQRTF (GRAV/WAVEN) 

OMEGAE=WAVEN*ABSF (CW+V) 

CALL ADMAB 

CALL COEFF 

CALL EXCITE 

FTWO--TI (2) 

EMIWO--TI(5) 

CALL MOTION 

L=(N/2)+1 

DELTA=DELTA*57.295779 
EPSIL=EPSIL*57.295779 

TNULL=TNULL/ (WA*WAVEN) 

ZNULL=ZNULL/WA 

IF(MINKRI) 16,16,15 

CALL BENDSH 

FORMAT (72H 

1 STATION) 

PRINT 921 

PRINT 923, KRIT 
SHNULL=SHNULL*BPL/ (GAMMA* (BSTAR (L) *WA) **2) 
BMNULL=BMNULL/ (GAMMA*WA*BSTAR(L) *BPL*BPL) 
ALPHA=ALPHA4 57 .295779 

BETA-BETA*57.295779 

DWL=WL/BPL 

FROUDE=V /SQRTF (GRAV*BPL) 

CONTINUE 

PRINT 903, FROUDE,DWL,OMEGAE, ZNULL, TNULL, BMNULL 
IF(MINKRI) 19,19,17 

IF(KRIT-MAXKRI) 18,19,19 

KRIT=KRIT+INCRES 

GO TO 15 


το) 





19 
20 


2 
22 


25 


IF(WL-BWL) 20,21,21 

WL=WL+DELWL 

TO TO 12 

IF (ABSF(V)-ABSF(VMAX)) 22,25,25 
=V+DELV 

GO TO 11 

CONTINUE 

GO TO 1013 
END(1,0,0,0,0,0,0,0,0,1,0,0,0,0,0,) 
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SUBROUTINE COEFF 04/18 


SUBROUTINE COEFF 

DIMENSION Y(21),TR(6),11 (6) ,ADD141[14),B EN (4),C55.(4),UK([ -)pUL(G), 
1DMASS1(21)»,QUANT(21,30),SKLAM(21)»BSiA (21)»C0XFSI1(21),SXf15T( +), 
2CTFST(2i1),STFST(21),XICLO 10) , DIXCZ210, N?XI(21)5,DA5FTCZC 223, 0OWEIGH(21) , 
AJSECOE(Z1),ABAR(21,30) ;WL(3 ),V(.0),-R QL(21 30) ,ABA<N( 2 ,»3u0)1 
4ADMASN(21,30) 

COMMON Y4,J4, SYMPS,DXIs,;ADDA,-"CEB,»CGGC, Z-EAL,Z1MAG, TREAL »sTIMAG,ZNULL, 
1)TNULL,DELTA,EPSIL,TR,TI4,V$sU0MASS,QUA' ; T; OMEGAC, SKLAM, KRITF,RO, 58RAV, 
2BSTAR,CXFST,SXFST,ALPHA, SHXYULL,XI, 5L T.,,BMNULL; GAMMA, DIXK, M; WÀA,WAVEN 
3,CW,ENOX1,SIGMA, TAU, FNULL,"MNULL ,DR«F:'  DWHEIGH, SE COE, TMASS,N, UR, UI ; 
4ABAR,PI; SHREAL; SHIMAG,BMREAL,BMIMAG, Y IERT,BPL,FAEQL,; ABARN, 
SADMASN,NSTA,K 

SMALL A = ADDA(1) 

J=M 

DO 10 I=1,M 

Y(1)=QUANT(I yK) 

CALL SIMPS 

ADMAS=SYMPS 

ADDA(1)=ADMAS+TMASS 

CAPITAL A = ADDA(4) 

DO 21 I=1,M 

Y(1)=QUANT(I1yK)#(XI(I)*w*2) 

CALL SIMPS 

ADDA(4)=SYMPS+YNERT/GRAV 

SMALL D = ADDA(3) = CAPITAL D = ADDal ) 

DO 30 I=1,M 

Y(1)=QUANT(1,K)+*xI1(1) 

CALL SIMPS 

ADDA(2)=SYMPS 

ADDA(3)=ADDA(2) 

SMALL В - ВЕЕВ(1) 

DO 40 I=1,M 

ENOXI(1)=(GAMMA*GRAV+(ABAR[ I,K)u..2))/(OMEGAE+*>) 

Y(I1)=ENOXI(I1) 

CALL SIMPS 

ВЕЕВ(1)-5ҮМР5 

CAPITAL B = BEEB(4) 

DO 50 I=1,M 

Y(I)=ENOXI(I)*(XI(I)*w*#2) 

CALL SIMPS 

TEMPORARILY 

BEEB(4)=SYMPS 

DIX(1)=-QUANT(2,K)/(2.*DXI1) 

DIX(M)ZQUANT(N,K)/(2. *DXI) 

DO 55 I=2,N 

DIX(I)=(QUANT(I-1+K)—-QUANT(I+1 yK))/(2 . w|)XI) 

DO 56 I=1,M 

Y(1)=DIX(1)+xI1(1)+x1(1) 

CALL SIMPS 

BEEB(4)=BEEB14)-ABSFIV)*(2.*ADDA[ 2)» SY MPS) 

CAPITAL E = BEEB(3) 

DO 60 I=1,M 

Y(1I)=DIX(I1)*w*XI(I) 

CALL SIMPS 

ETHRE=SYMPS#(ABSF(V)) 

DO 63 I=1,M 


ET 





SUBROUTINE COEFF C4/18 


63 Y(I)=ENOXI(I)*XI(I) 
CALL SIMPS 
EDNE=SYMPS 
ЕТНО-2.% (АВ85Ғ(У) ) «АОМА5 
BEEB(3)=EQNE-ETHRE 


C SMALL E = BEEB(2) 
BEEB(2)=BEEB(3)-ETWO 
C SMALL C =CGGC(i) 


O0 7O I21,M 
TO Y(I)s8STAR(I) 


CALL SIMPS 
CGGC(1) 2 GAMMA«SYMPS 
E CAPITAL C = CGGC(4) 


DO 71 [=1,М 
71 Y(1)=BSTAR(I)+*(XI(1)»«2) 


CALL SIMPS 
CGGC(4)=GAMMARSYMPS-(ABSFI(V))eBLEB( ) 
E CAPITAL G = CGGC(3) 


DO 80 I=1,M 

80 Y(I)=BSTAR(I)=XI(I) 
CALL SIMPS 
CGGC(3)2 SYMPS«*«GAMMA 

E SMALL G = CGGC(2) 

CGGC (2) =CGGC(3)-(CABSF(V) )#HEEB(i) 
RETURN 
END(1,1,0,0,0,0,1,15,0,1,0, ^. ,0,0,0) 


JOB TIME = 0.37 MIN. 


299 





SUBROUTINE EXCITE U4/1 


SUBROUTINE EXCITE 

DIMENSION Y(21),TR(6),T1(6),ADDA(4),B E8(4),CG65.(4) ;UR( D5,UI (6); 
1DMASS(21),QUANT(21,30), SKL3M(21),8Si A (21), CXFSI C21), SXF STOUC10; 
ФСТЕЗТ (21), 5ТЕЗТ (21), Х1(21),О1Х(21), -МХ1(21),О<АЕТ 21), ПИЕТОН (21), 
3SECOE(21)»ABAR(21,30),W1L (3), V(20),FR QL(21,30) ,ABARN(2 20), 
4ADMASN(21,30) 

COMMON Y,J,SYMPS,DXI,ADDA,?EEB+CGGU,Z 'EAL,ZIMA5)TREAL,TIMAG,ZNULL » 
1TNULL,DELTA, EPSIL, TR, TI, V DMASS,QUA iT, OMEGA ,SK1 AM,KRIT,KO,GRAV, 
2BSTAR,CXFST, SXFST,ALPHA, SHNULL, XI,03: F'4,8MNULL, GAMMA, DIX, M, WV ,WAVEN 
353 CW, ENOXI, SIGMA, TAU; FNULL,zMNULL ;D«^Fi ,DWEIGH,S.:-COE, IMASS,N, UR, UI » 
ABAR,PI,;SHREAL,SHIMAG,BMREAL,BMIMAG,;, Y ЕКТ, ВРЕ, Е EQL,ÄBARN, 
SADMASN,NSTA,K 

J=M 

DO 90 I=1,M 

FKLAM=(GAMMA+BSTAR( 1 )-( (WNAVEN+*CW)** )=QUANT(T,K)).WA 

SKLAM(I)=ENOXI(1)-(DIX(I)*ABSF(V)) 

CXFST(I)=(FKLAM*SINF(WAVEN*XI(I))+(WNAVYENFCHNA)=SKLAM(I) 
1*COSFIWAVEN*XI(1)))*(EXPF(-WAVEN*DRAF:(1)*SECOS (1))) 

SXFST(I)=(FKLAM*COSF(NAVEN*XI(TI))-(sAVENFCASWA)=ESKLAMI(I) 
1*SINF(WNAVEN+*XI(1)))* (EXPF(-NAVEN*DRF7(I)*S-CO-(1))) 

90 Y(I)=CXFST(I) 

CALL SIMPS 

C FONE 
TR(2)=SYMPS 
DO 92 I=1,M 
92 Y(I)=SXFST(I) 
CALL SIMPS 
@ -F TWO 
TI(2)2-SYMPS 
DO 93 I=1,M 
93 Y(I)=CXFST(I)+*XI(1) 
CALL SIMPS 
C EMONE 
TR(5)=SYMPS 
DO 94 I=1,M 
94 Y(I)=SXFST(I)+*XI(I) 
CALL SIMPS 
е -EMTWO 

TI(5)=-SYMPS 

SIGMA=ATANF(-TI(2)/TR(2)) 

TAU=ATANF(-TI(5)/TR(S)) 

FNULL=SQRTF(TR(2)+**.2+TI1(2)xx* 2) 

EMNULL=SQRTFITI(5)**24+TR(5)..2) 

RETURN 

END(1,1,0,0,0,0,1,1,0,1,0, 0.0.0) 


JOB TIME = 01.20 MIN. 


-90- 
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101 
102 


103 
104 
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110 
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n TIME = 
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SUBROUTINE MOTION 

DIMENSION Y(21),TR(6),11 (6), ADUVA(4),B 1b614),C063.(4)s+UR( ),UI1(5), 
1DMASS1(21)»QUANT(21,30),SKLaM(21),BSi4 Г 1), ХЕЪТ (21), БХУТО +), 
2СТЕ5Т(21),5ТЕ5Т(21), ХІ(21),01Х(.41), N.XI(21)yD<=`FT(/ )yWcCcIGH(Z : ) > 
ЗЗЕСОЕ (21), АВАК (21, 30), и! (3 ),V(=20), FR QLÍ(21,30) ,APAÁRNIZ ,30), 
4ADMASN(21,30) 

COMMON Y,J,SYMPS,DXI;ADDA,-:ULEB,CGGC; Z- E5L,21MAU, IREAL, T IMAG, ZNULL , 
iTNULL,DELTA,EPSIL,TR,TI4,V$,"MASS,QUA «T, OMEGA , SKLAM, &KKIT ,ROS, 5RAV, 
2BSTAR,CXFST,SXFST,ALPHA, SHNULL;XI,9^ T..,8MNULL $45 :MNA,.UIX, M. NY, hnAVEN 
3; CW, ENOXI,;SIGMA, TAU,FNULL,C:MNULL , DR F:,DWEIGH$,5-COE, IMASS,N, UR, UI ;, 
4ABAR, PI; SHREAL,SHIMAG4,BMREAL,BMIMAG, Y .ERT,BPL,F 4EQL , ABA" N, 
5ADMASN,NSTA 

DO 105 J=1,4 

GO TO (100,101,.02,103), J 

I=1 

GO TO 104 

1+3 

GO TO 104 

1=% 

GO TO 104 

I=6 

TR(I)=CGGC(J)-ADDA(J)*(OMEJAE+=2) 

TI(I)=BEEB(J)+OMEGAE 

DO 110 1-2,3 

DO 110 K=1,2 

IF(K#I-4) 108,110,108 

IPK=I+K 

UR(IPK)=TR(I)*TR(K+3)-TI(I)*TI(K+3)- T<«(K)eTR(I+-)+11(K)<11(1+43) 
UI(IPK)=TR(I)e*T1I(K+3)+TI (1) TR(K+43)-T-«(kK)TI(]4 )-T1(K)=1R(1+3) 
CONTINUE 

DO 111 I=1,6 

UR(I)=UR(1)/10000000000. 

UICI)-UI(I)/10000000000. 

DENUM=UR (4) ##24U1(4) ##2 

ZREAL=(UR(5) #UR(4)4U1I(5) #U1(4))/DENUM 

ZIMAG=(UR(5) #UI(4)-UI(5) «UR( 4) )/ DENUM 
TREAL=(UR(3) #UR(4)4UI (3) #Ui(4))/0ENUM 
TIMAG=(UR(3)*UI(4)-UI(3)+U<(4))/DENUM 

ZNULLZ SQRTF(ZREAL **2*Z2IMAGs *2) 

TNULLZSQRTF(TREAL *«2*TIMAG* ж 2 ) 

DELTA=ATANF(ZIMAG/ZREAL) 

EPSIL=ATANF(TIMAG/TREAL ) 

RETURN 

END(1,1,0,0,0,0,151,051,0,0,0,0,0) 


0.83 MIN. 
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SUBROUTINE BENOSH 02/02 


SUBROUTINE BENDSH 

DIMENSION У(21),ТК (6) „Т1 (6), АОПОПА ( 4), ВЕЕВ | 4), СССС: 4), ИВ (6) ЗУТ (6), 
10МА55 (21) „ОЦАМТ 21), ЗКГАМ( 21) ) „ВТАК ( 21) ,СХЕЗТ( 21), 5ХЕЗТ (21), 
2СТЕЗТ4 21) „5ТЕЗТ (21), Х1 (21) ,О1Х (21), ЕМОХТ (21) „ОКАЕТ (21) ΟΗΒΙΘΗ(ΣΖΙ 
35ЕСОЕ4 21), АВАН (21), АСКЕА! ( 21) ΑΟΙΜΑΟί 21) 

COMMON YoJeSYMPS »D0XI,ADDA, BEEB,CGGC,ZREAL»ZIMAG,TREAL,TIMAG, ZNULL, 
ITNULL CELTA EPSILo TRTI V 0OMASS, QUANT,OMEGAE, SKLAM,KRITs,RO&GRAV, 
2BSTAR,CXFST,SXFST,ALPHA, SHNULL, XI,BETA,BMNULL,GAMMA,DIX, ML WA WAVEN 
3; CW,ENCXI, SIGMA, TAU, FNULL, EMNULL 5, DRAFT4,0WEIGH, SECOE, TMASSSN, UR,UI , 
4ABAR,PI,SHREAL,SHIMAG,BMREAL,BMIMAG, YNERT,BPL 

Х1(-1. ) **KRIT 

IF(X) 111,111,112 

JzKRIT 

GO TO 113 

JsKRIT*1 

JENSsMRIT*1 

DO 114 Iz1,JENS 

ACREAL(I)Ss(ZREAL *TREALSXI(CI) )e8 COMEGAE* *2) 

ACIMAG(I)z(ZIMAG*TIMAG*eXICI) )* COMEGAE* *2) 

CTFESTCI) sACREALC I) #QUANT(TI) 
1+ ( QUANT(I)*2. 9 ABSF(V) *0OMEGAEI)STIMAG-SKLAMEI)S (CCZIMAG*TIMAGS* 
2XI (I) ) «*OMEGAE) -ABSF(V) «TREAL)-(GAMMAeBSTAR(CID))I* (ZREAL*TRESLe XI (I )) 
3e0XFSTtI) 

ЅТЕЅТ(І) «АСІМАС( І) #ОЧАМТ (І) 
1-(QUAAT(1[)+2.+ABSP(V)+OMEGAE)*+TREAL+SKLAM(I )e( ( (ZREAL+TREAL+ 
2X1(1))"»0OMEGAE)+ABSF(V)*TIMAG)-(GAMMARBSTAR(1)) (721 МАС+ ТЕМЯСе ХЕ (Г) ) 
3eSXFST (tI) 

Y(I) s0TFST(I) 

CALL :SIMPS 

SHREAL :Ü. 

DO 200 I=1,J 

SHREAL=5HREAL+ACREAL (1) DMASS( 1) 

SHREAL=SHREAL~(X #ACREAL( JENS) #OMASS(JENS)/2.) 

SHREAL=SYMPS-(X*CTESTIKRIT+1)eD0X1/2.)+SHREAL 

DO 118 1=1,Ј 

Y([)=STFST(I) 

CALL SIMPS 

SHIMAG=D. 

00 250 1є1,Ј 

SHIMAG=SHIMÁG+ACIMAG( 1)*+DMASS( 1) 

SHIMAG=SHIMAG-X* AC IMAG(JENS)*DMASS( JENS) 420. 

SHI MÁAG3SYMPS-(X«STESTIOKRIT*1)90XI/2.)* SHLMAG 

ALPHAzATANF (SHIMAG/SHREAL) 

SHNULL=SORTFISHREAL«=2+SHIMAG= «= 2) 

00 122 1=1,Ј 

Y(1)3CTFST(I)#@(XI(I)—-XI(KRIT+1)) 

CALL SIMPS 

BMREAL=0O, 

00 300 I=1,J 

BMREAL=BMREAL+ACREAL(I)=OMASS(I)*#(XI(I)—-XI(JENS)) 

BMREAL=BMREAL*ACREAL (JENS) “ОМА55(4ЕМ5)“0Х1І/8. 

BMREAL=SYMPS+(CTF$ST(KRIT+1)*#DXI*DXI/8.)+BMREAL 

DO 126 I=1,J 

Y(I)=STFST(I)*(XI(I)-XI(KRIT+1)) 

CALL SIMPS 

BMI MAG=0. 
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SUBROUTINE BENDSH _ 02/02 12 


DO 350 I=1,J 
350 ВМ| МАСО-ВР! МАС + АС ТМАС( Г) “ОМА55ЕТ) е (Х1 (Г) -Х1 ( ДЕМ) ) 
BMI MAG=BMIMAG+#AC IMAG( JENS) #OMASS( SENS) #DOXI/8. 
BMI MAG=SYMPS+(STFST(KRIT+1L) *#OXI#OXI/8.)+BMIMAG 
129 BETASATANF( BMIMAG/ BMREAL ) 
BMNULL2SQRTF(BMREAL#®#2+BMIMAG##2 ) 
RETURN 
ЕМО(1,0,0,0.0,0,0,0,0,1,0.0,0,0,0) 


JOB ТІМЕ - .05 MIN. 


SUBROUTINE SIMPS 04/18 


SUBROUTINE SIMPS 

DIMENSION Y(21),TR(6),T11 (6), ADDA(4),B.-EB(4),C6G62 14),UK(* ),UI(O6), 
1DMASS(21) QUANT(21,+30)y ySKLAM(21) 5yBSiA“ (21) | u XFS[(21) SXFSI(¿ `) 5 
ФСТЕЗТ (21), 5ТЕЗТ (21), Х1( 21), О1Х(21), М.Х1 (21) ,Ох ЕТ (е 1), СИЕТОН(2,), 
ЗЅЕСОЕ (21), АВАК (21,30) , ні (3С), У (20) ,ҒЕ О (21, 30) , АВАВМ( 2, 30), 
Ф+АОМАЅМ (21,30) 

COMMON Y; J4SYMPS,DXI,ADDA,PEEB,CGGC,Z*:EAL$,21MAG, IREAL, TIMAG; ZNULL , 
ITNULL DELTA EPSILs;TR, TI;V4 DMASS, QUANT, OMEGAC, SKLAM,KRIT ,RO,GRAV; 
2BSTAR,CXFST,SXFST,ALPHA, SHNULL,XI4,BET4,DBMNULL,GYMMA, DIX, M; W3,WAVEN 
3, CW, ENOXIs3SIGMA, TAU; FNULL; . MNULL, DRYF ; ,DWEIGH4, SCLCOE, IMAS S, N, UR,UI y 
4ABAR,PI,SHREAL,SHIMAG,BMREAL ,BMIMAG, YIERT,BPL,FREQL,ABAKNG 
БАОМАЗМ, МУТА , К 

SUM=Y(1)+Y(J) 

KARL=J-1 

DO 151 I=2,KARL 

X=(-1.)**I 

IF(X) 151,150,150 

150 YtI)22.*Y(I) 
151 SUM=SUM+2.*Y(1) 

SYMPS=(DX1/3.)*SUM 

RETURN 

END{1,1,0,0,0,0, 1,1,0,1,0,!,0,0,0) 


JOB TIME 01.56 MIN. 


-93- 





APPENDIX D 


MODEL RADIUS OF GYRATION 


The procedure used to determine the radius of gyration of the models is 
described. No information was available on the full scale radius of gyration 
for the catamaran or ALVIN. Therefore the radius of gyration was assumed to 
be 0.25 of the length. Because of the small weight of ALVIN compared to weight 
of instrumentation, which had to be positioned at the center of gravity, the 
largest radius of gyration that could be obtained was 0.232 of the length. 

The catamaran model was constructed with built-in lead ballast, and even with 


a drastic alteration of the model, the largest radius of gyration was 0.189 of 


the length. 
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Equipment: 


2 springs of same stiffness k 


1 stopwatch 


string, tape, set of ballast weights 


Procedure: 


a) 


b) 


ο) 


d) 


e) 


f) 


g) 


Suspend the two springs from the clamps which are attached to convenient 
points on the ceiling. The distance between the two springs should be 
2L. 

Make sure the ship model, without ballast, rests on level keel (or at 
desired trim) when put in water, and determine the center of gravity 

of the model. Assume that the center of gravity coincides with the 
center of flotation. Now put all ballast in model symmetrically about 
Ehe C.G. 

Measure off the distance L on each side of the center of gravity of the 
ship model. L can be of Ns convenient length. 

Suspend the model from the springs, using string wrapped around the 
model and secured in place by tape as shown in the figure above. 

Force the model to an up and down motion (heave) and count with the 
stopwatch the number of oscillations per minute. Thus, determine fa. 
Force the model to a pitching motion and count with the stopwatch the 
number of oscillations per minute. Thus determine es It is not 
necessary to achieve a pure pitching motion, completely uncoupled from 
heave, because the period of heave does not affect the period of pitch. 
Let W be the weight of the ship and w the weight of each of the springs. 


The moment of inertia of the system can be approximated by: 
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ο οι κ του I UNA (13) 
S s s g 


where κ. and КЕ are the radius of gyration of the system and the model, 


respectively. 


Therefore: k - 1 k - L (14) 


The radius of gyration of the system is given by: 


f 
= H ; - 2k - 2k _L 
k. = t L, since f = M & a M k. (15) 


where M is the mass of the system and k the stiffness of each spring. 


Hence, 


2 


f 
ΙΙ; as 
P 


For a given weight of model, the heaving period fa remains constant. 
£ | 
So, we can change the E ratio by changing the pitching period. To 


f 
P 


increase P dispose ballast weights closer to C.G., making sure they 


are still symmetrically disposed about the C.G. To decrease at dispose 
f 
ballast weights further away from the C.G. When E equals the value 


f 
P 


given from the above relation, the ballasting procedure of the ship is 


over. 
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APPENDIX E 


SUMMARY OF EXPERIMENTAL DATA 


The following data is a summary of the experimental data obtained at the 
M.I.T. Ship Model Towing Tank. It represents data read from Sanborn Oscillo- 
graph tapes. 

The catamaran single hull and catamaran data are for both directly ahead 
and astern seas. (Tables V-VIII). Table IX is for ALVIN in directly ahead 
seas and Table X for directly astern seas. Table XI is ALVIN data in directly 
ahead seas, while in the recovery position. Table XII is ALVIN pitching data 


and the catamaran heaving data, while in the recovery position. 


νος. 





TABLE V 


Catamaran Single Hull Experimental Data, h. = 0.375 inches 





z (in) | 9 (deg) 





TABLE VI 


Catamaran Experimental Data, ho = 0.375 inches 





A(ft.) z (in) O (deg) z /h 9 (rad) 
ο ο ο ο = en 
— | — |-5-- Kh 
ο 
3.0 0.0400 0.067 
4.0 0.0250 х 0.284 
4.5 0.0375 : 0.301 
6.0 0.0875 0.333 


10.0 0.3000 0.775 





Zoos 





TABLE VII 





Catamaran Single Hull Experimental Data, h_ = 0.75 inches 


9 on 
A(ft.) z (in) 9, (deg) 2 / 








TABLE VIII 


Catamaran Experimental Data, ho = 0.75 inches 


O, (tad) 
9 (deg) z hg KE 


0.0334 








3.0 0.10 ° 0.1333 
4.0 0.05 . 0.0667 0.1780 
4.5 0.0625 ° 0.0834 0.2500 


6.0 0.200 ° 0.2670 0.6000 


0.8910 





КРО Ци оборота то. 0.8000 


τ 





TABLE IX 


ALVIN Experimental Data, Directly Ahead Seas, ho z 0.375 inches 


0 (rad) 
A(ft.) z (in) o, (deg) z fh. X 


0.0045 





0.1068 


0.4000 


1.2000 





4.8000 





ALVIN Experimental Data, Directly Astern Seas, ho = 0.375 inches 


0.08 
0.1245 
0.375 


2.12 


— indicates data not obtained 


-100- 


Е 


е енен 





TABLE XI 

Recovery Experimental Data Directly Ahead Seas, h. = 0.375 inches 
9 (rad) 

Kh, 

0.005 

0.178 

1.175 

0.600 

0.890 


4.45 





TABLE XII 





Recovery Experimental Data Directly Ahead Seas, h_ = 0.375 inches 


ALVIN 
9 (deg.) z fh. 











0.454 


0.356 


0.320 


0.853 





2.140 
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